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We use the concept of resonant tunneling to calculate the thermopower of molecular nanosys- 
tems. It turns out that the sign of the thermovoltage under resonant tunneling conditions depends 
sensitively on the participating molecular orbital, and one finds a sign change when the transport 
channel switches from the highest occupied molecular orbital to the lowest unoccupied molecular 
orbital. Comparing our results to recent experimental data obtained for a BDT molecule contacted 
with an STM tip, we observe good agreement. 
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I. INTRODUCTION 

Studies of the origin of a voltage or current in nanosys- 
tems in the presence of a temperature gradient are an ex- 
tremely interesting and promising area in the field of nan- 
otechnologies [ll-Q- There are several important possi- 
ble future applications in several areas of devices, among 
them the development of nanothermosensors (see for ex- 
ample Q), which is especially urgent for a number of 
technological processes and for research in biology con- 
cerning the functioning of life. 

However, different from the classical description of 
thermoelectric phenomena, which is already challeng- 
ing enough, the necessity to apply strictly quantum- 
mechanical methods in the realm of nano-objects makes 
the whole problem an extremely difficult one, and a 
proper theory for studying transport phenomena in the 
most general setup does not yet exist. However, for the 
description of most experimental realizations of thermo- 
electric transport through nanon-structures, one can for- 
tunately make some simplifying assumptions. Usually, 
one can consider the system to consist of two metal- 
lic structures, which are typically very good conductors 
and which we will call leads, that are spatially separated. 
Hence, there will be no current fiowing between the leads. 
Placing an active element like a molecule between these 
leads will thus induce a transport path and, when volt- 
age or temperature differences between the leads are im- 
posed, to thermoelectric phenomena Q . The coupling of 
the molecule and the leads will be of tunneling type, i.e. 
one can usually assume that this coupling is rather weak. 

Although this setup seems rather straightforward, its 
actual experimental realization is by no means trivial. In 



particular, a good control over if a molecule is attached 
to both leads at all and how good the coupling to the 
leads is, is very hard to obtain Furthermore, when 
using two extended leads on a common substrate, the in- 
troduction of temperature gradients poses a further chal- 
lenge. A certain breakthrough along these lines has been 
achieved recently by using an STM tip to controlled pick 
molecules on a metal surface and generate well-defined 
break junctions. Here, certain control over the coupling 
between the molecule and the leads can be performed by 
moving the tip, and a temperature gradient can easily 
be applied by heating the metal while keeping the tip a 
fixed temperature 

In the present paper, we develop a non-equilibrium de- 
scription of the stationary thermotransport through such 
a nanostructure, using the idea of resonant tunneling. 
The model is rather simple, but can be solved analyt- 
ically, and the theoretical results can be directly com- 
pared with experimental data. We find a rather good 
agreement, which we believe is in favor of this simple 
and straightforward model. 



II. RESONANT TUNNELING WITH APPLIED 
TEMPERATURE GRADIENT 



Resonant electron tunneling of particles through a sys- 
tem of double potential barriers is very sensitive to a po- 
sition of electronic states in a quantum well [Q]. This 
circumstance can be used for the effective control of the 
tunneling process. As a model of a double-barrier tun- 
neling system, we assume the structure with the energy 
profile shown schematically in Fig. [1] j7l-[ll|. The Hamil- 
tonian describing the tunneling of electrons through such 
a structure can be chosen in the form 
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FIG. 1: Schematic model to study thermoelectric phenomena for a molecule attached to leads. A shows the case where the 
LUMO provides the resonant level, and B corresponds to a HOMO resonant level. The energies of orbitals and bias AV^ are 
chosen to provide a setup for thermopower measurmcnt, i.e. no electrical current flowing. 



The first term of this Hamiltonian is 

k(T k(7 

It describes electrons in the left electrode (source) and in 
the right (drain). Because we are not interested in the 
detailed properties of the leads, we assume that these 
charge carriers can be taken to be quasi-particles, and 
ai,L,A"'k,L,a) respectively aI,K,^(%_^ ,,) are the creation 
(annihilation) operators for these quasi-particle in source 
respectively drain. The dispersions are in the same spirit 
given by SL/nik) = ^^fc^/2m£,/fl., where mi^/n denote 
the effective masses for the left and right lead. We will 
assume in the following ttil — m/j = m for simplicity. 

The Hamiltonian Hw describes the electronic states in 
the nano-object. It can be written in the form 

Hw = J2^'- "i"" + ' (3) 

a 

where a labels the single-particle levels of the molecule, 
and Hi denotes possible interactions. The energy in the 
well depends on the applied bias AT^ and can be writ- 
ten as Ea = fa — eo 7 AV, where is the bare energy 
of the resonant state in the quantum well, AV the po- 
tential drop across the molecule, eo > the elementary 
charge and 7 a factor depending on the profile of the po- 
tential barriers (for identical barriers, 7 — 0.5). Finally, 
the Hamiltonian Ht describing the tunneling of electrons 
through the barriers has the conventional form 

Ht^ E (^-(fc)«M,.«a+h-c.) (4) 

ka,5=L,R 

Here, Ta„{k) is the matrix element of tunneling from 
source respectively drain to and from the molecule. As 



usual, we assume that both are the same and that they 
also do not depend on the applied bias. 

When we apply a constant external bias across the 
system, a nonequilibrium steady-state electron distribu- 
tion will result. We assume that the electron distribution 
functions in the electrodes (source, drain) are equilibrium 
ones, i.e. Fermi functions, due to the large volumes of 
these reservoirs, but their chemical potentials and tem- 
peratures can be different. The chemical potentials usu- 
ally encode a voltage drop across the nano-region, hence 
in our model /i^ = /i -I- A/i, = /i and A/i — cqAV . 

The simple setup of the system in Fig.[T]makes the eval- 
uation of nonequilibrium properties comparatively sim- 
ple. The important quantity entering all formula is the 
density of states (DOS) yo(i?) for the local level in the 
presence of the leads [12|- To calculate it, we need the 
retarded Green's function Ga.a{E) [isj, from which we 
can obtain the DOS as 

p{E) = VlmG„,„(^), 

TT ^ — ' 

a 

The electron distribution function g{E) in the quan- 
tum well is essentially nonequilibrium. It can be deter- 
mined from the condition of equality of the tunneling 
currents through the source and the drain. Then the 
distribution function has the form 

9{E) = ^^[TL{E)fL{E)+rn{E)ME)] (5) 
T{E) = Tl{E)+Tr{E) , (6) 

where F^ and Tn are the tunneling rates for source (L) 
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and drain (R), given by the expressions [l^l 

T^iE) = Y.\^U''^\'SiE~eUk)) (7) 

TRiE) = ^|r„,(fc)|'(5(ii;-£fl(fc)) , (8) 

k(7 

and /l(-E') and fniE) are the quasi-particle distribution 
functions in the source and the drain, respectively. They 
have Fermi-Dirac form and read temperature difference 



fL/RiE) 



E - ^^L/R , 



where ks is the Bohzmann constant, and Tj^/r, are tem- 
peratures in the source and the drain, respectively. 
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FIG. 2: Comparison of experimental data (blue circles 0]) 
with theoretical curve (bvll2[) is demonstrated. 



III. DOUBLE BARRIERS 
THERMOSTRUCTURES FOR RESONANT 
TUNNELING 

With the above formula for the distribution function, 
one can straightforwardly evaluate physical quantitites. 
For example, the occupancy of the molecule can be de- 
termined with the help of expression (3] 

Ua^-- J dE g[E) ImG'(a, a, E). 

Moreover, the net current Jsd between source and drain 
through the molecule is given by the equation [T^] 



Jsd — 



I G{E) {fL- fR)p{E)dE, 



(9) 



where G{E) = T l{E)T r{E)/T{E). Since we assume that 
the contacts are of tunnel type, the transition rates F^, 
Fi are exponentially dependent on the barrier widths 
and heights. Correlation effects between the electrons 
in a nano-structure encoded in Hi can be taken into ac- 
count by means of p{E), too and will in general dra- 
matically modify the properties [IB] ■ The theoretical 
description of this situation is at present possible in the 
linear response regime only (see Mravile and Ramsak [l^ 
for an overview). 

At this point one can only proceed analytically by as- 
suming a low transparency of the barriers, i.e. T{E) — > 0, 
and a simple energy level structure, for example a quasi- 
particle description of the electronic structure based 
on first-principle methods. The density of states p{E) 
then only depends on the energy structure of the nano- 
structure via 



where 5{x) is Dirac's function and Ca a weighting factor 
for the energy level Ea- For the current we then obtain 
the formula 

sd = -T- ^ Ca G{Ea) [fhiEa) ~ fR{Ea)) ■ 



J. 



(10) 



The distribution functions fL/R.{Ea) are exponentially 
dependent on the energy E^- Thus, when \Ea — Ea'\ ^ 
kB T, where Ea' denotes a neighboring molecular or- 
bital, there will be one particular energy Ea for which 
\Ea — is minimal. With respect to this orbital, trans- 
port through all other will be exponentially suppressed, 
hence we need to keep only this particular orbital in the 
calculations. In this approximation, [10] can be reduced 
to the form G{Ea)[fL{Ea) ~ fR{Ea)\ — 0, respectively, 
for small but finite Ti^/r{E), 



fdEa)^ fR{Ea)^Q. 



(11) 



According to our definition, Ea is one of the energy levels 
of a BDT molecule which has the smallest distance from 
the chemical potential, which will either be the highest 
occupied molecular orbital (HOMO) or the lowest unoc- 
cupied molecular orbital (LUMO). That molecule level 
is, as noted before, shifted by the voltage AF by a value 
Ai5 = — eo7Ay — — 7A/^. For an asymmetric barriere 
we have 7 = TR{Ea)/T{Ea)- Note that since the tunnel- 
ing matrix elements for source and drain were chosen to 
be the same, this expression remains well-defined in the 
limit V]^ir{E) 0. The solution of [TT] then becomes 



eoAl/ 



(Ea^fi) AT 



1 



T r 



(12) 



where as final step we also assumed that F2^/^(i?) « 
^L/R for energies close to /u. 

Our [T^l can now be compared with experimental data. 
In particular, such a comparison with the experiments 
on BDT with STM [2] resuhs in AV < 0. Acording 
to our setup in Fig. [1] this means that Ea < fi- Thus, 
the energy level Ea seen in experiment is a HOMO, i.e. 
the condition shown in Fig.[Tj3. In fact, the conductivity 
of the tunneling structure is determined by holes. The 
detailed comparison of the thermo- voltage dependence on 
AV is demonstrated in Fig. [H In the case of AT < T, 
[13 can be approximated by a linear relation 



eoAV = (Ea - p)AT/T. 



(13) 
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FIG. 3: Dependence of AV on the distance between electrodes 
at different temperatures. Results of calculation are compared 
with experiment 01 . 



Note that the total distance between source and drain is 
d = ttL + Gfi + aBDT: where cedt is the diameter of the 
molecule BDT. 

This theoretical length dependence of the thermoelec- 
trical voltage can be compared with experimental data 
0. The result is shown in Fig. [3l In our calculations, 
we fixed the parameters by obtaining the best fit for 
AT = 20K, resulting in ~ 0.2 nm, cedt = 0.5 nm, 
7 = 5 nm-\ Bq ^ 2 ■ lO^^^ - fj,)/T = 9.54 fiY/K. 
The remaining curves for the other temperatures were 
then calculated using these model parameters. The gen- 
eral agreement between theory and experiment is very 
good, only at larger AT we observe stronger deviations. 
Since our model is comparatively simple and naturally 
contains various assumptions, such a deviation is not sur- 
prising but to be expected. It can be attributed to var- 
ious sources, ranging from neglected barrier height fluc- 
tuations to influence of additional molecular orbitals or 
even molecular vibrations [16nl8i] . 



From the comparison with the experiment, we get /x — 
Ea ~ 2meV. 

To study the dependence of thermoelectric effects on 
the distance d between substrate surface (source) and 
STM tip (drain), let us consider our double-barrier sys- 
tem as a simple resistor network. In this case. [TT] and [TUl 
can be reduced to 

^ C„ G{E^)[h{E^) - fn{E^)] = J^a, (14) 
n 

with an externally imposed current Jsd = Ay/i?o, where 
i?o is the external impedance of the whole measurement 
system. We note that if the conditions eo|Ay| <^ ksT 
and \Ea — <^ ksT are valid, [T3]can be rewritten as 



[E^ - Ai)AT/T - eoAV = pAV, 



(15) 



where p = 4:hkBT/eoRoCaG{Ea). In the case of an 
open circuit Rq = co, we have p = 0, and [T51 and [T51 are 
the same. The quantities F^ and Tn entering G{E) are 
exponentially dependent on the barrier widths, i.e. we 
can approximate them as 



^LM. = To exp(-7aL,_R), 



(16) 



where a^^ji are the widths of left and right barriers, re- 
spectively, and 7 depends on the barrier height. The bar- 
rier lengths can in the experimental setup be controlled 
by moving the STM tip [31 • Finally, relation [T5] can be 
written in the form 



cqAV 



[Eg, ~ /i) Ar 

l + B{d) T ' 



(17) 



where 



B{d) = So[exp(7aL) -I- exp(7afl)], 



Bn 



IV. SUMMARY 

A theoretical description of transport through nano- 
structures makes a full quantum-mechanical description 
of the system mandatory. In contrast to bulk materi- 
als, one cannot even adopt some semi-classical approach 
based on e.g. Boltzmann equations here. Since one also 
needs to take into account the inherent non-equilibrium 
situation in many cases, solving this problem has become 
one of the most challenging tasks in modern condensed 
matter theory. A certain simplification arises when one 
can use the concept of resonant tunneling. This is usually 
possible in weakly contacted nano-objects like molecules, 
and allows to quite accurately describe the thermoelec- 
tric phenomena in these systems. 

We have presented here the calculation of thermotrans- 
port through a BEDT molecule contacted with a metallic 
substrate and a STM tip via generation of a break junc- 
tion 0. In the limit of only weakly transparent barriers 
we were able to obtain an explicit formula for the volt- 
age drop across the molecule as function of temperature 
difference between substrate and tip. We found that the 
experimental data are well described by a resonant tun- 
neling process involving the HOMO of the molecule. 

Modelling the dependence of tunneling rates between 
the molecule and substrate/tip with a simple exponential 
ansatz, we were furthermore able to give a closed expres- 
sion for the dependence of the thermoelectric effect on 
the distance between tip and substrate. The compari- 
son to experiment could be done by extracting the rele- 
vant model parameters from one set of data for a fixed 
AT — 20K only, reproducing the remaining curves with 
good accuracy. Thus, our formula describes very accu- 
rately the transport through such a nano-object, pro- 
vided the coupling to the leads is sufficiently weak. 

There are, of course, several severe simplifications in 
the model. The possibly most relevant is the neglect of 
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interaction effects on the molecule, as well as molecu- 
lar vibrations and also vibrations of the whole molecule 
between the contacts. These additional features can be 
taken into account [iGllTsj . and under very simplifying as- 
sumptions also solved analytically [18|] . We therefore be- 
lieve that such results, may they seem simple or straight- 
forward, are nevertheless very important steps to enhance 
our knowledge about transport through nano-structures 
and can actually also serve as benchmarks to test more 
elaborate theoretial tools. 
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